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摘要 射电 辐射 机 制 , 尤其 是 射电 暴发 现象 的 相干 辐射 机 制 , 是 天 体 物理 中 最 复杂 、 争 议 最 多 的 电磁 辐射 机 制 . 
于 受到 多 重 物理 因素 相互 牵连 的 复杂 影响 , 相干 射电 辐射 机 制 的 理论 研究 存在 很 大 的 难度 , 长 期 以 来 在 等 离 
- 子 体 辐射 和 电子 回旋 脉 泽 辐射 这 两 类 相干 辐射 机 制 间 争议 一 直 不 断 . 近年 来 , 人 们 开始 尝试 将 粒子 数值 模拟 方 
> 法 应 用 于 相干 射电 辐射 机 制 的 研究 , 并 已 经 取得 了 一 些 积极 的 进展 . 本 文 将 着 重 介绍 近年 来 的 粒子 模拟 研究 工 
作 及 其 取得 的 主要 进展 , 并 对 现存 的 一 些 问题 和 困难 进行 简要 评述 . 

关键 词 等 离子 体 , 波 , 辐射 机 制 : 非 热 , 方法 : 数值 
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1 引言 的 射电 相干 辐射 机 制 主要 有 两 类 , 即 基于 电子 等 离 
射电 天 文学 始 于 20 世 纪 30 年 代 美 国 年 轻 的 子 体 频率 附近 朗 缪 尔 波 的 波 - 波 非 线 性 不 稳定 性 激 
信 工程 师 Jansky 的 开创 性 工作 . 二 战 后 该 领域 得 到 发 过 程 的 等 离 了 体 辐射 (plasma emission) 机 制 “ 和 
了 迅速 发 展 , 并 在 短 短 的 几 年 时 间 内 探测 到 了 各 基于 电子 回旋 频率 附近 横 电 磁 波 的 线性 电子 回旋 
各 样 的 宇宙 射电 源 , 像 太阳 、 银 河 系 中心 、 超 新 星 脉 泽 不 稳定 性 激发 的 电子 回旋 脉 泽 辐射 (Electron 
遗迹 和 射电 星系 等 . 特别 是 , 有 关 太 阳 射 电 暴发 现 Cyclotron Maser Emission, ECME) 机 制 苞 站. 
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象 的 观测 研究 发 现 , 它们 不 仅 具 有 高 度 的 非 热 特征 ， 等 离子 体 辐射 机 制 由 Ginzburg 等 由 首先 提出 ， 
即 其 辐射 亮 温 度 76 远 大 于 辐射 源 电子 的 动力 学 ! 热 图 1 的 流程 图 展示 了 其 基本 物理 图 像 . 其 中 至 少 包 


< 温度 元 , 而 且 经 常 呈现 具有 漂移 特征 的 动力 学 谱 和 括 4 个 基本 物理 过 程 回 : 首先 由 高 能 电子 束 流 直接 

w 基 波 - 谐 波 对 的 双 带 结构 . 这 意味 着 这 类 暴发 辐射 BR BABAR YR (L, 见 文献 [38] 的 3.2 节 ); 然后 激发 的 

现象 应 该 主要 由 相干 的 等 离子 体 集体 辐射 过程 产 朗 缪 尔 波 与 离子 声波 (S) 的 波 - 波 非 线性 相互 作用 

生 , 它们 可 以 通过 在 等 离子 体 特征 频率 及 其 谐 频 附 产生 频率 在 等 离子 体 频 率 wue 附 近 的 基 频 射电 辆 

近 的 等 离子 体 不 稳定 性 激发 上 3. 目前 被 广泛 研究 射 (T) (BIL +S > T); 同时 , BRR SB 
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处 “动力 学 ?一 词 , 来 自 英文 单词 “kinetic” 的 翻译 , 其 英文 原意 表示 “运动 引起 的 "、“ 与 运动 有 关 的 ”等 意思 , 应 用 在 物理 学 领域 ,主要 描 
述 “ 由 于 粒子 运动 引起 系统 时 间 演 化 ”等 相关 的 物理 现象 . 近年 来 , 也 有 一 些 作 者 用 “ 动 理 学 ”代替 这 一 传统 翻译 ， 以 便 与 “<dynamic” 一 词 的 翻译 “ 动 
力学 ” 相 区 分 , 不 过 在 中 文 表述 里 ,“ 理 ”通常 是 表述 各 种 “静态 "关系 的 词 , 如 “纹理 ”、“ 道 理 "、“ 原 理 ” 等 . 所 以， 在 本 文中 , 为 了 强调 “时 间 演 化 ”这 
一 “动态 ”特性 , 我 们 依然 采用 “动力 学 ”一 词 来 表述 “kinetic” 的 含义 . 


m 


。 AN Ir HE 
AEH 
v i 4, 


(^ HAMM 
VC TI nar IV = 


64 45 X x 学 报 33 


KIER PER A E APE c RUP 8 BIS] B ER I BEL’) 流 的 卫星 探测 结果 , 对 等 离子 体 辐射 机 制 的 适用 
(BIL +S > L/); 而 两 个 反 向 传播 的 朗 缪 尔 波 通 过 性 进行 了 量 的 分 析 检 验 工作 ， 发 现行 星际 III 型 


非 线性 波 - 波 耦 合 进一步 在 等 离子 体 二 次 谐 频 附 近 ”射电 暴 与 电子 束 之 间 确实 存在 较 好 的 相关 性 , 但 
产生 谐 频 射电 辐射 (Ta) (ENL + L’ 一 Ty). 是 关联 的 朗 缪 尔 波 强 度 却 远 低 于 等 离子 体 辐 冉 

通过 以 上 描述 的 基本 物理 图 像 不 难 发 现 , 在 。 ”机 制 所 要 求 的 非 线 性 激发 阔 值 , 同时 也 没 能 观测 
人 磁化 等 离子 体 中 驱动 等 离子 体 辐射 机 制 的 ü | 到 电子 束 平行 动量 分 布 具有 “是 够 强 ” 的 布 居 数 反 
能 主要 来 自 能 够 有 效 激 发 朗 缪 尔 波 的 高 能 电子 束 、” HERHIEP 10. 尽管 在 随后 的 90 年 代 ，Cairns 等 人 和 


流 沿 磁场 方向 上 的 平行 能 量 , 因此 要 求 高 能 电子 。 Robinson 等 人 发 展 了 一 系列 更 精致 化 的 等 离子 体 
束 流 有 具有 布 居 数 反 转 的 平行 动量 分 布 函数 ful)， ”辐射 机 制 , 提出 了 朗 缪 尔 波 的 临界 激发 过 程 和 背 
即 wjdf(wj)/duy > 0 (其 中 由 为 电子 动量 在 平行 景 密度 涨 落 导致 行星 际 III 型 射电 暴 的 随机 生长 理 
于 磁场 方向 的 分 量 ). 自 上 世纪 70 年 代 以 来 , 针对 1600719), 但 等 离子 体 辐射 机 制 的 基本 困难 并 没有 
行星 际 II 型 射电 暴 的 卫星 观测 , 结合 高 能 电子 束 ”得 到 根本 性 的 解决 四. 
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图 1 ASR MR Ee iit eR 基于 文献 [3] 的 图 2. 
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Fig.1 Flow diagram for the theory of plasma emission. Note that this diagram was built up based on the Fig. 2 in Ref. [3]. 


相 比 于 等 离子 体 辐射 机 制 , ECME 机 制 的 物理 背景 等 离子 体 环境 必须 是 强 磁 化 的 , 以 致电 子 回 

图 像 相 对 比较 简单 , 其 中 的 基本 物理 过 程 是 通过 ” ” 旋 频 率 高 于 横 电 磁 波 穿 过 等 离子 体 传播 的 截止 频 
线性 电子 回旋 脉 泽 不 稳定 性 直接 放大 、 激 发 回旋 率 , 即 当地 等 离子 体 频率 ( 亦 即 wee > Wye) S71; 另 
频率 wee 及 其 谐 频 附近 的 横 电 磁 波 辐射 a4 2. 而 且 ， 一 个 条 件 是 高 能 电子 束 具 有 垂直 动量 (wu ) 布 居 数 
POMERIA ENI 电 辐射 的 偏振 特性 既 依 赖 于 高 。 反 转 的 速度 分 布 (f(wj)), Bld f (u,)/du, > 0, 如 环 
能 电子 束 的 分 布 状态 , 也 与 等 离子 体 磁 化 环境 的 Jyp, TEPSI TRACY, SBP 
物理 参量 有 密切 联系 - 因此 既 可 以 呈现 为 X 模 主导 式 的 电子 速度 分 布 . 因为 , 高 能 电子 束 回 旋 共 振 条 
的 右 旋 偏振 , 也 可 以 是 O 模 主导 的 左旋 偏振 [517. 件 的 相对 论 效应 导致 驱动 回旋 脉 泽 不 稳定 性 的 自 
不 过 , 有 效 地 产生 ECMEB 需 要 两 个 前 提 条 件 : 一 是 由 能 主要 来 自 高 能 电子 在 垂直 于 背景 磁场 方向 上 
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动力 学 运动 方程 进行 直接 自 洽 数值 计算 的 方法 , 也 


迄今 为 止 , BCME 机 制作 为 地 球 极光 千 米 波 辐 
射 n4 2 和 木星 、 土 星 等 行星 射电 辐射 29 以 及 一 些 
系 外 行星 射电 辐射 27 天 等 行星 磁 层 射电 辐射 现象 
的 基本 辐射 机 制 已 经 得 到 了 广泛 验证 和 普遍 认可 ， 
取得 了 巨大 成 功 [ 和 19. 但 是 , 在 应 用 到 太阳 射电 暴 
发 现象 时 , 局 限于 早期 人 们 对 太阳 大 气 结 构 的 认 知 ， 


AEDES y EEH], 直到 本 世纪 初 , 基于 对 
园 非 均匀 结构 的 新 认 知 ， 


是 研究 空间 无 碰撞 等 离子 体 微观 动力 学 过 程 的 最 
好 工具 , 可 以 比较 完整 地 再 现 等 离子 体 的 波 - 
粒 相 互 作 用 、 波 - 波 相 互 作用 等 非 线性 物理 过 
fel 44) PIC 数值 模拟 的 代码 最 先 由 Bunemanl5l 
和 Dawsonl46 在 50 年 代 末 60 年 代 初 开发 , 用 于 研究 
等 离子 体 中 粒子 之 间 的 电磁 相互 作用 . 因此 , 设 
计 PIC 数 值 模 拟 的 初衷 就 是 为 了 模拟 等 离子 体 
中 的 动力 学 过 程 . 迄今 为 止 PIC 数值 模拟 仍然 


Wu 等 人 重新 检验 了 
ECME 机 制 应 用 于 解释 太 


和 行星 际 III 型 射电 暴 
的 可 能 性 B39. 近年 来 , Wu 等 人 提出 普遍 存在 于 
各 种 空间 等 离子 体 环境 的 满 动 阿尔 文 波 能 够 对 
ECME 机 制 产生 重要 影响 , 从 而 定性 地 改变 其 物理 
特性 [39, 推动 了 ECME 机 制 在 各 类 太阳 射电 暴发 现 
象 的 研究 中 的 应 用 B37. 在 此 基础 上 , Chen 等 人 及 
Wu 等 人 也 进一步 发 展 了 一 种 包含 高 能 电子 束 自 4 
阿尔 文 波 效应 的 新 型 ECME 驱 动机 制 E7, 35-391. f 
够 比较 好 地 克服 传统 ECME 机 制 应 用 到 各 类 太 虹 
射电 暴发 现象 时 所 遇 到 的 主要 困难 (详情 可 参见 本 
专辑 中 唐 建 飞 等 人 的 综述 文章 mo). 

不 过 , 一 方面 由 于 相干 辐射 机 制 的 工作 效率 比 
较 敏 感 地 依赖 于 辐射 源 区 的 等 离子 体 环境 , 而 我 
们 目前 对 太阳 磁化 等 离子 体 大 气 物 理 参 量 及 其 分 
布 状态 的 了 解 尚 远 远 达 不 到 足以 判别 不 同 辐射 机 


LH oe [tt 


制 所 需要 的 精准 程度 ; 男 一 方面 , 由 于 相干 辐射 机 
制 通常 比较 容易 受到 其 他 未 知 波 模 等 复杂 因素 的 


干扰 和 影响 , 给 理论 模型 的 研究 也 带 来 了 很 大 的 
不 确定 性 . 因此 , 近年 来 的 一 些 研究 开始 引入 对 模 
型 简化 过 程 依赖 性 较 低 的 粒子 模拟 方法 , 即 所 谓 
的 Particle-in-Cell (PIC) 模 拟 , 对 相干 射电 辐射 机 


是 模拟 再 现 无 碰撞 等 离子 体 中 各 种 动力 学 物 
理 过 程 最 方便 的 计算 工具 之 一 . PIC 数值 模拟 
的 基本 理论 和 应 用 在 经 典 参 考 教科 书 中 有 很 好 
的 解释 ， 如 Birdsal 等 上 i、Hockney 等 上 以 及 其 他 
有 关 PIC 数 值 模拟 的 综述 文章 5 53. Birdsall 在 文 
献 [53] 中 对 PIC 数 值 模拟 代码 的 发 展 过 程 进行 了 回 
顾 . 

PIC 数值 模拟 的 核心 就 是 数值 求解 麦克 斯 
*B 77 fe 2A ((1)- (4) X) DAC Y B) 1 fe 2 d 2] 77 
程 ((5) 式 ): 


V.E(g,t) = 22" | (1) 
€0 

V- B(z,t) =0, (2) 

vx Ez.) =F). (3) 

V x B(z,t) = ia T(z, t) + Moeo DES ; (4) 

ED — De v) s) x BE] 6 

它们 描述 了 电场 及 PABA RF Xe BECBÉ I 


制 进 行 模拟 研究 , 并 取得 了 一 些 积极 的 进展 . 在 接 
下 来 的 第 2 节 中 , 我 们 将 对 PIC 模拟 方法 进行 简要 
介绍 . 随后 的 第 3、4 节 分 别 介绍 弱 磁 化 和 强 磁 化 等 


离子 体 环 境 中 相干 射电 辐射 机 制 的 PIC 模拟 研究 及 
其 主要 结果 . 第 5 节 进 一 步 对 比分 析 弱 磁化 和 强 磁 


化 等 离子 体 情形 下 PIC 模拟 研究 的 结果 , 并 在 最 后 
的 第 6 节 进 行 简 单 的 总 结 和 讨论 
2 PIC 数值 模拟 方法 简介 

所 谓 PIC 模拟 , 是 一 类 基于 等 离子 体 粒子 微观 


AE 
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间 t、 空 间 z 的 演化 规律 . 其 中 , gy、ma、p 以 及 J 分 
别 为 粒子 带电 量 、 粒 子 质 量 、 电 荷 密度 以 及 电 
Tit, 而 Lo 和 eo 则 分 别 为 真空 介 电 常量 和 真空 磁 导 率 . 
Vay 等 缴 总 结 了 PIC 数 值 模拟 中 所 涉及 的 各 种 数值 
计算 方法 . 

PIC 数值 模拟 方法 的 基本 思想 如 图 2 所 示 . 首 
先 , 模拟 区 域 被 划分 为 许多 网 格 单 元 , 电磁 场 的 
信息 仅 分 布 于 网 格 (单元 的 交叉 ) 点 , (E, B); (其 
中 7 代表 不 同 的 网 格 点 ),， 而 计算 粒子 (Ma) 则 可 处 
于 模拟 区 域 的 任何 位 置 ( 见 图 2 (a)). 然后 , 采用 “ 插 
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值 (interpolation)” 计 算 方法 来 获取 粒子 任意 位 置 
上 的 电磁 场 信息 , 并 通过 插值 得 到 的 电磁 场 来 进 
一 步 计 算 粒 子 在 其 位 置 上 所 受 的 洛 伦 效 力 瓦 (其 
中 i 代表 不 同 的 计算 粒子 ), 见 图 2 (b) 中 的 步骤 1 以 
及 在 该 洛 伦 兹 力作 用 下 , 粒子 于 新 时 刻 的 速度 元 和 
frr ( 即 图 2 (b) 中 的 步骤 2). 接 下 来 , 再 通过 对 
新 时 刻 粒 子 位 置 、 速 度 信息 (元 全 ;的 插值 , 可 得 到 


(b) 


Boris push 
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Integration of eqs of 

motion, moving 
particles 

Ue E E 


Weighting 
(E, B); zx E 
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该 时 刻 各 个 网 格 点 上 电荷 与 电流 密度 (o, 妃 ; 的 信 
息 ( 即 图 2 (b) 中 的 步骤 3), 再 将 该 电荷 与 电流 密度 
信息 代入 麦克 斯 韦 方 程 组 进行 数值 积分 , 以 获得 新 
时 刻 的 电磁 场 信息 (EE,B); ( 即 图 2 (b) 中 的 步 又 4). 
整个 PIC 数 值 模拟 即 是 对 上 述 4 个 步骤 的 循环 往复 ， 
直至 模拟 结束 . 


J assignment Weighting 


(x, v); =e (p, J)i 


Maxwell eqs 


Step 4 
Integration of field eqs 
on grid 
(E, B); - (p. J); 


Force interpolation 


图 2 


Fig.2 Panel (a): schematic diagram of a PIC simulation domain. Panel (b): main loop in a PIC codel 


3 ” 弱 磁 化 等 离子 体 中 相干 射电 辐射 

机 制 的 PIC 模 拟 研究 

对 相干 射电 辐射 机 制 、 尤 其 是 ECME 机 制 影 
响 最 显著 的 背景 等 离子 体 参量 是 电子 回旋 频率 
和 等 离子 体 频率 的 比 : Re = wee /woe. TERI HRe 
< 1 的 弱 磁 化 等 离子 体 环境 中 , ECME 机 制 将 因为 
ECME 的 激发 频率 (~ wee) 低 于 横 电磁 波 传播 的 截 
止 频 率 (w。) 而 不 能 有 效 工作 . 因此 , 相干 射电 辐射 


X|(a): PIC 数值 模拟 区 域 的 示意 图 . 


(b): PIC 数值 模拟 代码 的 基本 运作 流程 4. 
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粒子 动量 和 电磁 场 具 有 3 分 量 ) 的 PIC 模拟 方法 , JE 
究 了 电子 束 激发 二 次 谐 波 (2f,) 的 强度 对 模拟 维度 
以 及 对 不 同 物理 参量 的 依赖 性 . 图 3 展示 了 该 研究 
的 主要 结果 , 其 中 第 1 (2) 行 为 1.5 (2.5) 维 模拟 结 

图 3 (a)、(b) 列 则 通过 电场 分 量 轧 , ((a) 列 ) 和 磁 
场 分 量 B,((b) 列 ) 的 能 量 在 平行 于 背景 磁场 的 色散 
空间 (kj 一 w) 的 分 布 来 展示 等 离子 体 中 静电 波 (ES) 
和 电磁 波 (EM) 的 不 同 波 模 的 激发 情况 ( 注 : ,和 
Bz 分 别 平行 和 垂直 于 背景 磁场 , 其 主要 分 别 存在 


的 ECME 机 制 只 能 在 频率 比 Re > 1 的 强 磁化 等 离 
子 体 环境 中 有 效 地 产生 射电 其 发 辐射 , 而 等 离子 体 
辐射 机 制 原则 上 不 受 这 个 磁化 条 件 的 限制 . 这 里 ， 
我 们 先 介绍 弱 磁 化 等 离子 体 中 相干 射电 辐射 机 制 
的 PIC 模拟 研究 , 下 一 节 再 介绍 强 磁 化 等 离子 体 中 
相干 射电 辐射 机 制 PIC 模拟 研究 的 情况 . 

针对 弱 磁 化 等 离子 体 (wc。= 0.01wpe) 背 景 ， 


于 等 离子 体 静 电波 和 电磁 波 中 ), (c)-(f) 列 则 展示 了 
等 离子 体 中 几 个 被 较 强 激发 的 静电 波 和 电磁 波 (如 : 
BAZ ARYL. Be el Fe HG BB AR YL! SH) A He (RK 
发 强度 ) 对 背景 离子 质量 mi (Q7) ae Ba ih 
FET; (7) HAF RRR FEV, ((e 列 ) 以 及 柬 流 电 
子 数 密度 nb ((f) 列 ) 的 依赖 . 该 研究 发 现 静电 二 次 
谐 波 (ES 一 2f,, 图 3 (a) 列 ) 的 产生 不 依赖 于 模拟 空 


Kasaba 等 时 分 别 采 用 1.5 维 (一 维 空间 坐标 、 粒 子 
动量 和 电磁 场 具 有 3 分 量 ) 和 2.5 维 (两 维 空间 坐标 、 


间 的 维度 , 而 电磁 二 次 谐 波 (EM 一 2f) 仅 在 二 维 
的 模拟 中 被 激发 ( 见 图 3 (b) 列 ). mE, 反 向 传播 朗 
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间 等 离子 体 中 相 


ChinaXiv 合 作 期 刊 


F 射电 辐射 过 程 的 粒子 模拟 研究 3 期 


一 2f,) 的 激发 ( 见 图 3 (e). (£) 9D). 基于 二 次 谐 波 的 
激发 强度 对 不 同 物理 参数 的 依赖 性 以 及 满足 3 波 


(ES — LF, 可 代表 离子 声波 ) 的 功率 . 这 意味 着 反 向 
传播 朗 缪 尔 波 的 产生 可 能 是 正 向 传播 朗 绢 尔 波 经 
由 离子 声波 的 感应 衰变 过 程 ( 即 L + S 一 也) 产生 
的 , 而 不 是 经 由 朗 缪 尔 波 自发 衰变 过 程 ( 即 L 一 L 
十 S) 产 生 的 , 因为 后 者 一 般 不 依赖 于 离子 声波 的 功 
度 ( 访 ) 和 束 流 电 子 数 密度 (nb) 


率 . 另外 , 电子 束 流速 
的 增加 都 有 利于 正 、 反 向 朗 缪 尔 波 (L、L)、 二 次 
次 谐 频 电磁 波 (EM 


谐 频 静电 波 (ES 一 2f) 以 及 二 


图 3 电场 分 量 忆 。((a) 列 ) 和 磁场 分 量 B。((b) 列 ) 的 能 量 在 平行 于 背 


E 1 1 1 1 
eo 10000m, 1600m, 100m,.0004T, (1/3)T, T, 


-120 
co 10000m, 1600m, 100m, .0004T, (1/3)T, T, 3T, SV, 10v, 


( 波 1、 波 2、 a ee 
十 wa = ws, ky + ko = ka, Kw, Rl ki, i = 1,2,3 2) 
别 为 相互 作用 3 个 ; al ORAM a ces Kasa- 

ba P4 — 2p TR HH: 二 次 谐 频 静 电波 产生 于 两 个 正 
5 — EMBL +L ES — 2f,); 
fi eet EAS BE. Jc] BRR CHR 6r DU] A 
产生 二 次 谐 频 电磁 波 ( 即 L +L 2 EM 一 2f,) 的 基 
fili. 


m 


3T, SV, 


景 磁场 的 色散 空间 (ks 一 w) 的 分 布 . (c)-(£) FU: Mo RUBIA KUE(L). Jen] 
传播 朗 缪 尔 波 (L')、 频 率 小 于 0.1wse 的 低频 静电 噪音 (ES 一 LF)、 二 次 谐 频 静电 波 (ES 一 2fy) 以 及 二 次 谐 频 电磁 波 (EM — 2f, ) 的 振幅 对 背景 离子 
质量 mi ((c) 列 )、 背 景 离子 温度 T((d) 列 )、 电 子 束 流速 度 Vi((e) 列 ) 以 及 束 流 电子 数 密度 nl，((f) 列 ) 的 依赖 . EIT: 展示 1.5 维 模拟 结果 ， 下行: B 
示 2.5 维 模拟 结果 . 该 图 摘自 文献 [54]. 


Fig.3 The parallel (to the background magnetic field) dispersion (i.e., kz — w) diagram of the electric E; (column (a)) and 


magnetic B, (column (b)) components. Columns (c)-(f): amplitude of beam-excited Langmuir wave (L), backward propagating 


Langmuir wave (L’), electrostatic low-frequency noise below 0.lw,. (ES — LF), electrostatic secondary harmonic wave 


(ES — 2f,), and electromagnetic secondary harmonic wave (EM — 2f,) dependence on background ion mass m; (column (c)), ion 


temperature Ti (column (d)), electron beam velocity V, (column (e)), and electron beam number density np (column (f)). Top 


row: results from 1.5D simulations. Bottom row: results from 2.5D simulations . See Ref. [54] for the original plots. 


在 空间 射电 辐射 的 观测 中 , 除了 基 波 和 二 次 
谐 波 被 探测 到 以 外 , 与 之 相伴 的 更 高 次 谐 波 也 经 
常 被 捕捉 到 , WT, MER BAS FE = 
R55) 在 地 球 激 波 前 兆 区 中 甚至 出 现 过 五 次 谐 
波 [59. 通过 2.5 维 PIC 模 拟 , Rhee 等 [6 验证 了 漂移 


速度 达到 0.5c (c 为 真空 中 光速 ) 以 上 的 高 能 电子 束 
对 非 磁 化 等 离子 体 中 高 次 谐 频 波 的 激发 . 图 4 通 
xl nda EE, (图 (a)) 和 磁场 分 量 B。( 图 (b)) 的 能 
量 在 平行 于 背景 磁场 的 色散 空间 (ks — w) 的 分 布 
来 分 别 展示 所 研究 的 等 离子 体 中 静电 波 和 电磁 


chinaXiv:202307.00115v1 


64 4 大 


波 的 激发 情况 . 他 们 指出 : n 次 谐 频 电磁 波 (nH) 的 
形成 是 由 静电 朗 织 尔 波 (L) 与 (n 一 1) 次 谐 频 电磁 
波 ((n 一 1)H) 的 相互 看 合 过 程 , BL + (n 一 1)H 一 
nHl61-63， 而 n 次 谐 频 静电 波 (nL) 则 产生 于 静电 朗 
ZERIE(L)5(n- 1) 次 谐 频 静电 波 ((n 一 1)L) 的 相 


C k/ Ge 


1.2e-10 


2.4e-02 1.7e-10 
图 4 ”电场 分 量 。( 图 (a)) 和 磁场 分 量 B。( 图 (b)) 的 能 量 在 平行 于 背景 磁场 的 色散 空间 (kj — w) 的 分 布 . 该 图 摘 
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HA, BIL + (n — 1)L > mLI63. 近来 Yao 等 64 通 
过 2.5 维 PIC 模拟 研究 了 背景 等 离子 体 数 密 度 的 不 
均匀 性 对 高 能 电子 束 激 发 基 、 谐 频 电 磁 波 的 影响 ， 
结果 显示 背景 密度 的 不 均匀 性 也 有 利于 更 高 次 谱 
频 波 的 产生 . 


文献 [60]. 


Fig.4 The parallel (to the background magnetic field) dispersion (i.e., ks — w) diagram of the electric E, (panel (a)) and 


magnetic B; (panel (b)) components. See Ref. [60] for the original plots. 


受 早 期 计算 资源 缺乏 的 限制 ， Kasaba 等 5 
和 Rhee 等 6 的 PIC 模拟 时 长 大 约 只 有 200/wpe (对 
应 于 5/wpi, wpi 是 离子 等 离子 体 频率 ), 该 时 长 太 短 
无 法 探测 到 低频 离子 声波 的 激发 . 此 外 , 这 些 PIC 
模拟 所 涉及 的 空间 尺度 也 较 小 (512AXp。x 512Ape, 
Ape 为 电子 的 德 拜 长 度 ). 由 于 PIC 模拟 可 涉及 的 最 
小 波 数 反 比 于 模拟 空间 的 尺度 (pain x 1/ D), 进而 
模拟 中 较 小 的 空间 尺度 也 不 利于 对 具有 小 波 数 等 


降低 PIC 模拟 中 离子 的 数值 噪声 , 是 通过 PIC 模拟 
正确 研究 等 离子 体 辐射 机 制 理论 的 前 提 . 随 着 计算 
机 技术 不 断 的 发 展 以 及 计算 资源 的 不 断 优化 , 在 后 
来 的 相关 模拟 研究 中 , 模拟 时 长 、 模 拟 空间 尺度 、 
网 格 粒子 数 确实 都 得 到 了 有 效 的 提升 . 

Umedal6] 通 过 2.5 维 PIC 模拟 对 比 研 究 了 单 、 
(对 称 的 ) 双 向 电子 束 对 频率 比 Re = 0.1 的 弱 磁 化 等 
离子 体 中 二 次 谐 频 电磁 波 的 激发 以 及 不 同 背 景 离 


离子 体 波 进行 探测 分 析 . 同时 , 每 个 网 格 中 的 粒子 
数 也 较 少 (20-80 个 粒子 每 网 格 ). 而 PIC 模拟 中 每 个 
网 格 包 含 的 模拟 粒子 数 较 少时 会 增强 等 离子 体 中 
离子 的 数值 噪声 ( 非 物理 的 扰动 )@, 65-661, 这 些 数 
值 噪声 可 能 进一步 加 强 离子 对 朗 缪 尔 波 的 散射 以 
KERR Vel BA A BA EB E. 因此 , 较 短 的 
模拟 时 长 、 较 小 的 模拟 空间 尺度 以 及 较 少 的 网 格 
粒子 数 都 严重 干扰 相干 射电 辐射 机 制 的 PIC 模拟 研 
究 的 准确 性 . 因此 , 正如 Umedale] 指 出 的 , 有 效 地 
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子 温度 对 单 束 流 等 离子 体 中 二 次 谐 频 电磁 波 激 发 
的 影响 . 该 研究 的 各 模拟 中 所 采用 的 参数 及 其 区 
p| ILKI, HPna qi/ Mir Va; 40V. at 39] 73 55 BT 
体 中 3 种 不 同 粒子 ( = [1, 2,3]) 的 粒子 数 密度 占 比 、 
(被 正 电 子 荷 质 比 归 一 化 的 ) 荷 质 比 、( 被 电子 热 速 
度 归 一 化 的 ) 漂 移 速 度 和 (被 电子 热 速度 归 一 化 的 ) 
热 速度 . 模拟 A 和 模拟 B、 模 拟 C 分 别 为 包含 双 和 单 
向 电子 束 的 等 离子 体 , 具有 单 向 电子 束 的 模拟 B 与 
模拟 C 则 具有 不 同 背 景 离子 温度 . 图 5 则 展示 了 各 被 
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模拟 等 离子 体 中 电场 的 各 分 量 的 能 量 随时 间 的 演 子 束 情形 中 , Bc ty BE AR Be PY EE JS Wl fe Fe HOR 
化 , 其 中 电场 分 量 , ( 忆 ,) 平 行 (垂直 ) 于 背景 磁场 可 ”线性 不 稳定 性 直接 激发 , 其 能 量 要 高 于 单 向 电子 束 
反映 等 离子 体 中 静电 波 ( 电 磁 波 ) 的 特性 . 他 们 的 研 。” 情形 中 由 非 线 性 波 - 波 克 合 过 程 产生 的 反 向 朗 缪 尔 
FRM: 与 等 离子 体 辐 射 理论 预期 的 一 致 , (模拟 A UK, 这 也 与 等 离子 体 中 非 线 性 过 程 的 激发 效率 一 般 
中 由 对 称 的 双向 电子 束 所 激发 的 ) 能 量 相当 的 正 、 远 低 于 线性 不 稳定 性 激发 效率 的 理论 预期 一 致 . 此 
反 向 朗 缪 尔 波 的 波 - 波 相互 作用 可 有 效 增强 等 离子 Sh, 在 离子 温度 较 低 (全 < 区) 的 单 向 电子 束 等 离子 
体 中 电磁 波 的 激发 ( 见 图 5 (a) 中 电场 分 量 B 的 峰值 。” 体 中 (模拟 B), 反 向 朗 缪 尔 波 可 能 是 由 正 向 朗 缘 尔 
能 量 远 高 于 具有 单 向 电子 束 的 等 离子 体 即 图 5 (c))  ” 波 通 过 衰变 产生 的 , 即 L > L+S, 而 当 离 子 温度 与 
另外 , 单 向 电子 束 激发 电磁 波 的 增长 率 (7) 也 要 小 电子 温度 相当 时 (了 ~ T, EAC), 由 于 离子 热 扰动 
于 双向 电子 束 激发 电磁 波 的 增长 率 ( 见 图 5 (a) 中 的 的 增强 , 反 向 朗 缪 尔 波 则 可 能 主要 由 离子 声波 与 正 
wpet = 50 — 2501 图 5 (c) Fwpet = 200 — 4000: 7R Ve] BJZR ZR IR HEHE, 即 L+S 一 L ( 见 图 5 
的 电场 B: 的 非 线 性 增长 阶段 ). 这 一 结论 与 随 (b) 中 wpet = 400 — 900). 

后 Ganse 等 68 -709 的 研究 结果 一 致 . 因为 在 双向 电 


Ez 


Parametric decay 


0 100 200 300 400 500 600 700 800 io 200 400 600 800 1000 
to. to, 


Run B 
y ~ 0.005 for run B and 0.012 for run C 
107 
0 200 400 600 800 1000 
to. 


图 5 图 (a): BL ArnTBIS3A2HE(E,. ETE.) etc oa E BENE AE. 图 (b) 和 (c): BEB AUG C rh Fy op at E ELE HIRE s PEE B E] TRU PI TR 
化 . 该 图 摘自 文献 [67]. 


Fig.5 Panel (a): time history of energy densities for electric field components E,, E,, and E; in run A. Panels (b) and (c): time 


history of energy densities for electric field components E, and E, in run B and run C, respectively. See Ref. [67] for the 


original plots. 


表 1 不 同 模拟 中 的 参数 分 布 . 该 表 摘 自 文 献 [67]. 


Table 1 Parameters for different simulation runs. See Ref. [67] for the original table. 


Run mi api Va Va m gq/m2 Va V2 m q3/ma Vas Vis 
A 0.98 —1 0 1 0.01 —1 25 1 0.01 —1 —25 1 
0.99 —1 0 1 0.01 —1 25 1 1 0.000625 0 0.00625 
C 0.99 —1 0 1 0.01 —1 25 1 1 0.000625 0 0.025 
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在 相干 射电 辐射 理论 的 PIC 模拟 研究 中 , 降低 
数值 噪声 对 所 研究 波 模 的 影响 一 直 是 研究 者 需要 
重点 关注 的 问题 , 虽然 较 多 的 网 格 粒 子 数 可 以 有 效 
降低 数值 噪声 , 但 也 将 大 大 增加 计算 资源 的 负担 ， 
延长 完成 一 个 模拟 工作 所 用 的 计算 时 间 . 因此 , 采 
多 少 网 格 粒 子 数 最 为 适合 , 需要 进行 收敛 性 测 
YA, 即 对 比 采 用 不 同 网 格 粒子 数 的 模拟 结果 是 否 随 
粒子 数 增加 而 渐 趋 稳定 . Thurgood Xt EE dt fT 
了 收敛 性 检验 , 他 们 采用 2.5 维 PIC 模 拟 方法 研究 了 
网 格 粒子 数 为 125、250、500、1000 的 差别 , 结果 
发 现 粒 子 数 为 500 和 1000 时 计算 结果 的 差异 在 可 接 
受 的 范围 内 . 同时 , 他 们 的 计算 结果 也 显示 , 相对 
T ABR JI PRR, 离子 声波 增长 到 可 明显 
探测 和 证 认 的 水 平 需要 更 长 的 时 间 . 男 外 , 他 们 还 
注意 到 在 以 往 的 模拟 研究 中 , 大 都 采用 了 有 具有 较 
高 相对 密度 的 高 能 电子 束 , 如 n/no > 1% 的 情形 ， 
其 中 mo 是 背景 热电 子 的 数 密 度 . 而 在 大 多 数 空间 环 
境 ( 如 太阳 风 ) 中 , 实际 观测 到 的 高 能 电子 束 密 度 往 
往 更 低 . 有 鉴于 此 , 他 们 进一步 对 比 研究 了 低 密 度 
(nby/no = 0.5796, 模拟 1) 弱 电子 束 和 高 密度 (n/no 
= 596, 模拟 2) 强 电子 束 对 等 离子 体 中 电磁 波 的 激 
R. 表 2 展 示 了 各 模拟 所 采用 的 物理 参数 , 其 中 wb 
与 Aub = 公分 别 为 束 流 电子 的 漂移 速度 与 束 流 电 
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表 2 低 (模拟 1)、 高 (模拟 2) 密 度 高 能 电子 束 流 的 模拟 参数 . 
该 表 摘自 文献 [71]. 


Table 2 Simulation parameters for tenuous (run 
1) and dense (run 2) energetic electron beam. 
See Ref. [71] for the original table. 


Run  n»/no wp/Avb = vv/Ve Notes 
1 0.0057 16 Weak beam 
2 0.05 8 Dense beam 


在 上 述 的 各 项 PIC 模拟 研究 中 都 忽略 了 一 个 重 
要 问题 , 即 没有 考虑 在 初始 AURA TRE 
a 子 东 所 携带 的 束 流 电流 可 和 能 引起 背景 等 离子 体 
响应 而 被 “电流 中 性 ”化 的 效应 . 实际 上 , 高 能 电子 
束 的 “注入 ”有 一 个 逐渐 增强 的 过 程 , 随 着 束 流 电流 
的 增加 , 必然 导致 背景 等 离子 体 的 电磁 响应 而 激发 
形成 “返回 电流 ”, 以 致 有 效 地 抵消 束 流 电流 引起 的 
强 磁 扰动 效应 . Henri 等 中 在 他 们 的 2.5 维 PIC 模拟 
中 采用 密度 为 必 = 10n% 的 反问 慢 束 电子 成 份 来 中 
和 高 能 电子 束 引 入 的 束 流 电流 , 且 仅 考虑 mb = 
5x10-* — 5x10-3no 的 低 密 度 高 能 电子 束 对 等 
离子 体 中 电磁 波 的 激发 . 图 7 展示 了 该 研究 的 
主要 结果 , 其 中 图 7 (a)-(d) 展 示 了 被 不 同 的 束 
流 电 子 的 漂移 速度 (vw。， 图 (a) 和 (b)) 和 束 流 电子 
的 数 密度 (nb/no， 图 (c) 和 (qd)) 所 激发 的 正 向 朗 纪 
f 


子 的 热 速度 . 而 图 6 则 通过 平行 于 束 流 运动 方向 的 
电场 分 量 。((a) 和 (b) 列 ) 和 垂直 于 模拟 平面 的 磁 
场 分 量 B, ((c) 和 (d) 列 ) 在 平行 于 束 流 运动 方向 的 色 
散 空 间 (ks-w) 的 能 量 分 布 来 分 别 展 示 等 离子 体 中 
静电 波 和 电磁 波 在 不 同时 段 (第 1 行 : 50 < wpet < 
150, 第 2 行 : 450 < wpet < 550, 第 3 行 : 850 < wpet 
< 950) 的 激发 情况 (图 中 入 p 为 电子 德 拜 长 度 ). 结果 
显示 , 在 低 密度 电子 束 情形 中 被 激发 的 静电 波 主要 
T BEAR I E (LEG. (a) 列 ), 而 在 高 密度 电子 束 中 朗 
缘 尔 波 和 束 流 模 都 能 被 激发 , 但 束 流 模 占 主导 成 
份 , 被 激发 的 归 缪 尔 波 的 能 量 较 低 , 难以 通过 非 线 
性 波 - 波 耦 合 有 效 激 发 辐射 电磁 波 ( 见 图 6 (b) 列 ). Al 
此 低 密 度 电子 束 可 以 更 有 效 地 激发 电磁 波 辐射 ( 见 
图 6 (c) 和 (d) 列 ). 


REL Bc In] B] 2 AR JE (L^) WA — UR Ve A E R 
波 (Ta2p) 的 能 量 (Vw) 随 时 间 的 演化 . Re. Un 
别 为 束 流 电子 的 热 速 度 、 束 流 电 子 的 初始 动能 
fp = wp/(27)=wpe/(27) 为 电子 等 离子 体 频 率 . 图 7 
(e) 和 (f) 则 展示 了 二 次 谐 频 电 磁 波 的 能 量 在 不 同方 
向 (其 中 0° 代 表 束 流 电 子 运动 的 方向 ) 上 的 分 布 对 不 
司 的 束 流 电 子 数 密度 和 漂移 速度 的 依赖 . 图 7 (g) 则 
通过 垂直 于 模拟 平面 的 磁场 分 量 B, 的 能 量 在 色散 
空间 — w 的 分 布 来 展示 等 离子 体 中 电磁 波 的 激 
发 情况 . 这 些 模拟 结果 显示 : 不 同 于 高 密度 电子 束 
Ti 309. 低 密 度 束 流 电子 的 数 密度 和 漂移 速度 对 
二 次 谐 波 激发 效率 的 影响 不 大 ( 见 图 7 (aj-(d)), 但 
被 激发 二 次 谐 波 的 方向 性 则 强烈 依赖 于 束 流 电 子 
的 数 密度 和 漂移 速度 ( 见 图 7 (e) 和 (f)). 且 总 体 上 来 
Vt, 低 密度 高 能 电子 束 激发 的 电磁 波 以 二 次 谐 波 为 
主 ( 见 图 7 (g)). 
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6 (a) 和 (b) 列 : 模拟 1 ((a) 列 ) 和 模拟 2 ((b) 列 ) 所 获得 的 3 个 不 同时 间 段 的 电场 分 量 瓦 。 的 能 量 在 色散 空间 (KE。 — w) 的 分 布 . (c) 和 (d) 列 
与 (a) 和 (b) 列 相似 , 但 (c) 和 (d) 列 展示 了 磁场 分 量 刀 = 的 能 量 分 布 情况 . 该 图 摘自 文献 [71]. 


Fig.6 Columns (a) and (b): the ks — w dispersion diagram of the electric field component E; for run 1 (column (a)) and run 2 
(column (b)) over three time periods. Columns (c) and (d) are similar to Columns (a) and (b) but for magnetic component B+. 


See Ref. [71] for the original plots. 


在 磁化 等 离子 体 中 , 磁场 能 够 引起 不 可 忽略 的 。 ”GL)、 束 流 模 (Beam-Langmuir, BL)、 基 波 (O-F)、 
混杂 效应 . 在 非 磁化 等 离子 体 中 , 通常 只 有 静电 波谱 波 电磁 模 (也 J、 哨 声波 (W) 和 2Z 模 波 (Z-mode, LZ) 
和 横 电 磁 波 这 两 种 扰动 模式 存在 . 而 背景 磁场 的 引 ”的 激发 情况 (注意 : 不 同 列 中 波 数 k 和 频率 的 范围 
入 将 可 能 导致 大 量 混杂 扰动 模式 的 出 现 , 如 哨 声 。。 的 不 同 ). 该 研究 发 现 : 在 可 辨认 的 被 束 流 电子 所 激 
波 、 混 杂 波 等 这 些 混杂 滤 模 的 出 现 也 势必 对 磁化 。 发 的 波 模 中 不 仅 有 静电 的 朗 缪 尔 波 和 束 流 模 (GL 和 
过 离子 体 中 相干 射电 辐射 过 程 产生 重要 影响 , 尤其 。 Br. 见 图 8 (a) 列 ) 以 及 基 波 和 谐 波 电磁 模 (O-F 和 HH 
在 各 种 空间 和 天 体 等 离子 体 环境 中 , 背景 磁场 几乎 。 见 图 8 (b) 和 (4) 列 ) 而 且 还 有 磁化 等 离子 体 中 经 常 


是 无 处 不 在 的 . 近来 ，Chen 等 中 通过 2.5 维 PIC 模 ee 

是 无 处 不 在 的 , LOR, Chen 等 ”通过 2.5 维 PIC 樟 oy smi pe Rr MESE (WHILZ, 见 图 8 (c) 和 (b) 列 ) 
拟 详 细 研 究 了 高 能 电子 束 流 对 弱 磁 化 进一步 的 分 析 发 现 在 所 有 被 激发 的 波 模 中 HOR 
(Rr = 0.1) 等 离子 体 中 各 种 波 模 的 激发 情况 . 图 8 展 。 BOTA, 人 人 人 
示 了 不 同 的 电磁 场 分 量 (B; = E. Ej En Bes 模 的 能 量 占 主导 , 而 且 被 激发 的 朗 比尔 波 以 及 基 、 
B, 或 B,) 在 沿 O4B = arctan(kz/ky) = 30° (k, fk, 谐 频 电磁 波 的 能 量 与 束 流 模 的 能 量 存在 正 相 关 性 ， 
分 别 为 波 和 撩 在 xz、z 方 向 上 的 分 量 , x 方向 平行 于 这 意味 着 束 流 模 对 朗 缪 尔 波 以 及 基 、 谐 频 电磁 波 
背景 磁场 ,> 方向 垂直 于 模拟 平面 ) 的 色散 空间 的 的 产生 过 程 具 有 重要 贡献 , 甚至 后 者 诸 波 有 可 能 直 
fie Bt 4) f OK c H 7 YK (Generalized Langmuir, 接 产 生 于 束 流 模 的 衰变 过 程 ， 
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图 7 图 (a) 和 (b): 正 (L)、 反 向 传播 (L') 朗 缪 尔 波 以 及 二 次 谐 频 电磁 波 (T2t) 的 能 量 随 时 间 的 演化 . 图 (c) 和 (qd): 工 十 和 Ts2t 的 能 量 随 时 间 的 演化 . 


图 (e) 和 (f): Ts 在 不 同方 向 上 的 能 量 分 布 . 图 (g): 磁场 分 量 B, 的 能 量 在 色散 空间 天 


ww 的 分 布 . 该 图 摘 


文献 [72]. 


Fig.7 Panels (a) and (b): time evolution of the wave energy of the forward (L), backward (L’) propagating Langmuir waves, 


and harmonic electromagnetic emission (T;). Panels (c) and (d): time evolution of the energy of the L + L’ and Taof. Panels (e) 


and (f): energy distribution of T»; along different directions. 


See Ref. [72] for the original plots. 


Ni 等 47 引用 具有 垂直 动量 反 转 的 DGH 
(Dory-Guest-Harris) 分 布 5g 的 高 能 电子 , 1E 99 B 
化 (Rr = 0.1) 等 离子 体 环境 中 进行 了 2.5 维 的 PIC 模 
拟 研究 , 发 现 DGH 分 布 可 以 有 效 驱 动 弱 磁 化 等 
离子 体 中 的 双 等 离子 体 共 振 不 稳定 性 57-79， 引 
起 上 混杂 模 (Upper Hybrid，UH)、2Z 模 以 及 哨 声 
模 的 激发 . 这 些 首先 被 激发 的 波 模 之 间 的 波 - 波 
相互 作用 可 引发 其 、 谐 频 电 磁 波 辐射 . 图 9 通过 


Panel (g): the 天 一 ww diagram of the magnetic field component B+. 


频 电 磁 波 的 产生 依旧 是 基于 非 线性 波 - 波 相互 作用 
过 程 . 因此 , Ni 等 "4 的 PIC 模 拟 计算 模型 依然 可 
以 归 类 于 等 离子 体 辐射 机 制 的 范畴 . 

图 10 展 示 了 3 个 不 同 PIC 数 值 模拟 研究 [7 73-74 
中 所 获得 的 高 能 电子 束 流 激发 等 离子 体 中 不 同 波 
模 的 能 量 ( 被 高 能 电子 初始 能 量 归 一 化 ) 随 时 间 的 演 
化 . 结合 图 5、7、10 不 难 发 现 : 通过 非 线性 波 - 波 耦 
合 激发 的 基 、 谐 频 电 磁 波 的 能 量 不 仅 远 小 于 高 能 


Manley-Rowei KEMXA Aki + ko = 局 来 展示 
基 、 谐 频 电 磁 波 辐射 生成 的 可 能 的 原因 : 相反 方 
向 传播 的 Z 模 波 和 W 模 波 的 耦合 可 激发 基 频 电磁 波 
辐射 , 即 Z +W > OF ( 见 图 9 (c)), 而 相向 传 所 
的 UH 模 波 的 耦合 则 可 产生 谐 频 电 磁 波 辆 射 , BIUH 
+ UH’ > H ( 见 图 9 (d)). 与 传统 等 离子 体 辐射 机 
制 的 理论 框架 不 同 , 这 里 双 等 离子 体 共振 不 稳定 性 
的 驱动 自由 能 来 自 电子 的 垂直 自由 能 , 但 基 频 和 谐 


B T REIHE BERE, 也 远 小 于 其 他 静电 或 混杂 波 模 
的 能 量 . 意味 着 通过 非 线 性 波 - 波 厢 合 产生 辐射 电 
人 磁 波 具有 极 低 的 转化 效率 . 综合 上 述 的 PIC 模 拟 研 
完结 果 不 难 发 现 , 在 Re < 1 的 弱 磁 化 等 离子 体 环 境 
中 , 辐射 电磁 波 都 是 经 由 非 线性 波 - 波 相互 作用 而 
产生 的 , 尽管 参与 非 线 性 耦合 的 波 模 可 以 是 静电 波 
或 混杂 波 . 这 与 基于 等 离子 体 不 稳定 性 的 理论 预期 
也 是 一 致 的 , 因为 在 弱 磁 化 等 离子 体 环境 中 , 能 够 
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H 


通过 线性 不 稳定 性 直接 放大 的 电磁 波 频率 低 于 辐 TAR, 因而 会 通过 波 粒 相互 作用 将 获得 的 能 
射电 磁 波 传播 的 截止 频率 , 上 述 PIC 模 拟 的 结果 也 ”返还 给 粒子 , 从 而 降低 辐射 电磁 波 的 激发 效率 . 这 
确实 显示 没有 这 类 辐射 电磁 波 出 现 . 不 过 需要 强调 ”也 可 能 是 导致 最 终 的 辐射 电磁 波 能 量 过 低 的 原因 
的 是 : 在 上 述 所 有 的 PIC 数 值 模拟 中 都 采用 了 周期 ” 之 一 . 

性 边界 条 件 , 这 意味 着 电磁 波 被 束缚 于 等 离子 体内 


(a) BL and GL (b) LZ and O-F (OW (d) H 
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图 8 "ARE, (AUT). E,B217). E. (8317) EWO, 5 —arctan (k,/k,) = 30° 的 色散 空间 的 能 量 分 布 来 反映 束 流 - 朗 缪 尔 波 和 广义 朗 缪 尔 

波 (BL 和 GL,，(a) 列 )、2Z 模 波 和 基 频 电磁 波 (LZ 和 O-F,(b) 列 ) 的 激发 情况 , 而 磁场 分 量 B, (GALT). B (第 2 行 )、B;( 第 3 行 ) 用 于 展示 哨 声 波 (W/， 

(c) 列 ) 和 谐 频 电 磁 波 (H，(d) 列 ) 的 激发 情况 . 各 图 上 欠 加 展示 的 曲线 对 应 于 冷 磁化 等 离子 体 理论 中 的 不 同 波 模 . 在 (a) 列 的 各 图 中 , 竖 直 的 虚线 则 对 
应 于 束 流 模 的 色散 曲线 . 该 图 摘自 文献 [73]. 


Fig.8 Dispersion diagrams along direction 0,5 = arctan(k,/k,) = 30? of the electric field components E, (row 1), E, (row 2), 
E, (row 3) to present the excitation of the BL and GL (column (a)), LZ and O-F modes (column (b)), magnetic field 
components B, (row 1), By (row 2), B; (row 3) are presented for the excitation of the W (column (c)) and H modes (column 
(d)) of waves. Overplotted lines are corresponding dispersion curves of wave modes in magnetized code plasmas. The straight 


dashed line in column (a) represents the beam mode. See Ref. [73] for the original plots. 
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基 、 谐 频 电 磁 波 的 产生 可 能 经 历 的 波 - 波 相互 作用 过 程 . ERU EL SCA [74]. 


Fig.9 The possible configurations of wave-wave interaction processes for the generation of fundamental and harmonic 


electromagnetic waves. See Ref. [74] for the original plots. 


(a) , 0 (c) 
10 1(b) Case R 10-3 
10 a 10: E, + E, (BL) 1074 
E i 10 3 1075 
= 10+- 4 49-34 Ey + Ez (GL 
& oW 1071 | uos 
2 "uU | 10 E 1077 
o [neti 10-5 3 1078 
E | 
AT] w: j 10-5 j E+B(H) E (O-F) 10-9 
aa 10-7 E 10710 | 
10° 10-8 dH 4 
200 400 600 800 0 500 1000 1500 2000 0 1000 2000 3000 
ta. t (w52) t (wpe) 


] 


FBO P~ 1077179. 图 (b): 


电子 的 动能 演化 ( 黑 线 ). 电场 分 量 的 能 量 演 化 


X 10 不 同 PIC 数值 模拟 所 获得 的 高 能 电子 束 流 激 发 等 离子 体 中 不 同 波 模 的 能 量 (被 高 能 电子 初始 能 量 归 一 化 ) 随 时 间 的 演化 . 图 (a): 束 流 、 背 景 执 
红线 , MESTRE, Ey. Ez). 磁场 分 量 的 能 量 演化 ( 蓝 线 , ~ 10 ^ DIETS B, WB. AB, Jt 


ER Vt- BAB AVE (BL), 哨 声 波 (W), BEA (O-F) D REB REDE (H) US. 图 (c): 上 混杂 波 (UH), 哨 声波 (W)， 
Z 模 波 (Z), 基 频 (O-F) 以 及 谐 频 电磁 波 (H) 41. 


Fig.10 Temporal profiles of energies (normalized by the initial kinetic energy of energetic electrons, i.e., Exo) of various wave 
modes excited by energetic electron beam in different PIC simulations. Panel (a): the black lines show the kinetic energy 
associated with the beam and background electrons, the red lines are for the electric field energy (from x, y, and z components 
when reading from top to bottom) and blue lines present the magnetic field energy (B, component at ~ 1074, and B, and By, 


components at noise level ~ 1077) (741, Panel (b): BL for the beam-Langmuir mode, W for whistler, O-F for the fundamental, 


and H for the harmonic plasma emission"?), Panel (c): UH for upper hybrid mode, W for whistler mode, Z for Z mode, O-F for 


O mode around the fundamental plasma frequency, and H for harmonic plasma emission 


4 强 磁 化 等 离子 体 中 相干 射电 辐射 


机 制 的 PIC 模拟 研究 


在 频率 比 Re > 1 的 强 磁化 等 离子 体 环 境 中 , 电 


子 回旋 脉 译 不 稳定 性 能 够 直 
子 体 截止 频率 的 辐射 电磁 波 ， 


Œ 


接 放大 频率 高 于 等 离 


因而 相干 射电 辐射 能 


[74] 


够 由 线性 不 稳定 性 直接 放大 产生 , 这 就 是 相干 射电 
辐射 的 ECME 机 制 . 不 过 , 等 离子 体 辐射 机 制 中 入 
非 线性 波 - 波 耦 合 等 基本 物理 过 程 原则 上 在 强 磁化 
等 离子 体 环境 中 也 可 以 发 生 , 而 且 产 生 电 磁 波 的 频 
率 也 可 以 在 截止 频率 之 上 . 因此 , 强 磁 化 等 离子 体 
环境 中 相干 射电 辐射 机 制 的 PIC 模拟 研究 将 有 助 于 


cr 
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我 们 进一步 厘清 等 离子 体 辐射 机 制 和 ECME 机 种 
的 相对 重要 性 . 

强 磁化 等 离子 体 中 相干 射电 辐射 最 典型 的 例 
子 就 是 来 自 地 球 极 区 磁 层 的 极光 干 米 波 辐射 . 基 
于 地 球 极 区 不 同 高 度 处 的 等 离子 体 参 数 , Pritch- 
ett 等 人 采用 1.5 维 B99 和 2.5 维 I81-89 的 PIC 模 拟 方法 ， 
研究 高 能 电子 束 对 Re € [10, 40] 的 强 磁 化 等 离子 体 
中 电磁 波 的 激发 情况 . 他 们 的 模拟 研究 显示 , 在 
这 类 强 磁化 等 离子 体 环境 下 电子 回旋 脉 泽 不 稳定 
性 确实 能 有 效 地 激发 电子 回旋 频率 附近 横 电 磁 波 
的 增长 , 而 且 且 有 非常 高 的 激发 效率 , 可 以 将 大 
约 7% 的 电子 束 能 量 转 换 至 辐射 电磁 波 . 特别 是 , 被 


Ke 8 
激发 的 电磁 波 模 由 反常 X 模 所 主导 , 与 Wu-Lee 模 


im 


[sj 


FH 


周 晓 伟 等 : 空间 等 离子 体 中 相 了 


F 射电 辐 射 过 程 的 粒子 模拟 看 
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-二 
= 


间 和 天 体 等 离子 体 环 境 中 , 这 种 高 密度 电子 束 的 情 
形 还 是 比较 少见 的 , Pritchett 等 go -5 的 模拟 结果 
对 低 密 度 高 能 电子 束 的 情形 是 否 依然 成 立 , I IF 
进一步 的 研究 . 

取 高 能 电子 束 数 密度 为 nb/mno = 5%、 频 率 比 
Ret = 5，Lee 等 &3-89 利 用 2.5 维 PIC 模拟 研究 了 有 具有 
环 - 束 动量 分 布 的 高 能 电子 束 在 强 磁 化 等 离子 体 
中 对 高 频 电磁 波 的 激发 效率 . 他 们 的 结果 显示 : 被 
激发 的 主要 波 模 有 高 频 的 Z 模 、X 模 以 及 较 低 频 的 
fic BH (Beam-Langmuir, B-L) 和 束 流 - 哨 声 
波 (Beam-whistler, B-W), 但 激发 能 量 在 不 同 波 模 
之 间 的 分 配 比较 复杂 , 不 仅 依赖 于 激发 波 模 的 传播 


型 9 的 理论 预期 很 好 地 一 致 . 同时 , 他 们 也 对 比 
分 析 了 周期 性 (无 限 , infinite) 边 界 条 件 和 吸收 (有 
WR, bounded) 边 界 条 件 以 及 (考虑 到 极光 干 米 波 辐 


方向 和 环 - 束 电 子 的 投射 角 加 , 而 且 也 依赖 于 环 - 束 
电子 的 能 量 7 QE: 此 处 7 为 环 - 束 电 子 的 特征 能 量 ). 
图 12 展 示 了 有 具有 不 同 特征 能 量 和 投射 角 的 环 - 束 电 


射 源 区 位 于 极 区 磁 层 中 的 磁场 存在 不 均匀 性 这 
一 观测 事实 ) 磁 场 梯度 对 电子 回旋 脉 泽 辐射 激发 
效率 的 影响 . 表 3 和 图 11 清 楚 地 展示 了 这 些 模拟 与 
常规 电子 回旋 脉 泽 辐射 模拟 中 (通常 采用 周期 性 


子 对 等 离子 体 中 波 的 激发 的 区 别 , 其 中 (a)-(d) 列 
通过 不 同 电磁 场 分 量 在 沿 05 (= arctan(k,/k,), 
Ex. Jra FT HAT BS 84 39)7; [Hl 
E Rf kw II rn B Se E 4r dn OR HE as SE ALT Men 
有 哪些 波 得 到 激发 , (a) 7) WI zs T A-R E T f 


ae 


边界 条 件 以 及 均匀 密度 、 磁 场 ) 电 子 回旋 脉 泽 辐 


射 激发 效率 的 区 别 . 表 3 中 席 为 地 球 的 半径 、Q。 为 
电子 回旋 频率 、A 为 一 个 正方 形 模拟 网 格 的 边 长 
而 为 整个 (正方 形 ) 模 拟 区 域 的 边 长 . 图 11 (a) 通 
过 比较 周期 性 边界 条 件 和 吸收 边界 条 件 的 模拟 


(bounded) 边 界 条 件 的 模拟 中 非 热 电子 损失 更 多 
直方 向 上 的 自由 能 以 用 于 激发 电磁 波 辐射 . 图 11 
) 则 对 比 了 均匀 、 非 均匀 磁场 的 模拟 中 电磁 波 的 
量 (Paa) 随 时 间 的 演化 的 差别 , FP Pin AVE 
电子 的 初始 动能 . 结果 如 上 文 提 及 的 : 吸收 边界 
条 件 ( 见 表 3 以 及 图 11 (a)) 和 磁场 在 垂直 磁场 方向 
上 的 梯度 ( 见 图 11 (b)) 都 有 利于 提升 电子 回旋 脉 译 


E 
ER. 


的 变化 量 以 及 不 同 波 模 所 携带 的 能 量 随时 间 
的 演化 . 通过 图 12 我 们 可 以 发 现 : 当 电 子 
较 高 时 主要 激发 X 模 基 波 (X1) 的 电磁 波 辐 
12 第 1 行 ), 而 能 量 较 低 的 电子 束 更 有 利于 Z 模 基 波 
(2Z1) 的 激发 , O 模 基 频 波 (O1) 和 X 模 谐 波 (X2) 的 激 
发 非常 微弱 , X1 波 则 没有 被 激发 ( 见 图 12 第 2 行 ). 另 
外 , 他 们 也 注意 到 : 投射 角 越 大 的 环 - 束 电子 越 有 
利于 Z1、X2 以 及 O1 模 的 激发 . 上 述 结果 似乎 与 线 
性 电子 回旋 脉 泽 不 稳定 性 的 理论 预期 不 太一 致 . 
Lee 等 [83- 史 认为 该 结果 可 能 是 由 于 初始 增长 较 快 
的 静电 波 模 的 非 线 性 效应 导致 环 - 束 电子 自由 能 的 
重新 分 配 , 以致 呈现 出 比较 复杂 的 波 能 分 配 关系 . 


用 273 
不 过 , 在 他 们 的 模拟 中 并 没有 考虑 初始 时 刻 环 - 束 


重 


辐射 中 电磁 波 的 激发 效率 . 不 过 , 需要 指出 的 是 ， 
在 Pritchett 等 Bo -83 的 模拟 工作 中 高 能 电子 束 的 数 
密度 占 总 电子 数 密度 的 主要 成 份 , 高 达 总 电子 密度 
75% ( 见 表 3). 在 极光 电子 加 速 区 该 情形 可 能 是 
合理 的 , 因为 在 磁 层 区 域 源 自 磁 层 的 高 能 电子 数 通 
常 远 远 多 于 源 自 电 离 层 的 热电 子 . 但 是 , 在 其 他 空 
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电子 携带 电流 的 中 性 化 问题 , 尤其 是 投射 角 较 小 的 
环 - 束 电子 , 较 强 的 初始 束 电 流 可 引发 电流 不 稳定 
性 可 能 导致 静电 波 模 反常 的 快速 增强 , 而 抑制 了 其 
他 电磁 波 模 的 正常 增长 . 这 也 可 能 是 导致 其 模拟 结 
果 中 出 现 激发 自 


能 被 异常 地 重新 分 配 的 原因 之 
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表 3 无 限 边界 和 有 限 边 界 PIC 模 拟 中 电子 回旋 脉 泽 辐 射 的 效率 . 该 表 摘自 文献 [80]. 


Table 3 Efficiency of electron cyclotron maser instability in infinite and bounded PIC simulations. See 
Ref. [80] for the original table. 


Plasma Width 
Altitude RE wpe/Q。 Efficiency in Infinite System (96) Efficiency in Bounded System (96) 
Ls/A Lswpe/c 


Pure Magnetospheric Population 


1.5 0.050 4.5 32 tel 11.1 

64 2.3 14.6 

128 4.6 15.4 

256 9.1 13.7 

512 18.3 12.9 

768 27.4 11.9 

1024 36.6 11.2 
1.75 0.030 6.0 256 5.7 12.1 

alis 2.0 0.025 2.9 512 9.1 8.6 


Y— 2.5 0.050 0.5 256 9.1 2.8 


75% Magnetospheric, 25% Ionospheric 


1.75 0.030 1.5 256 5.7 5.0 


(b) 
Nonuniform B 
0.12 
Bounded 
> 3 Uniform B 
[qv] g 0.08 
— a 
c 
- o 
z 
= a 
0.04 
0 
0 1600 3200 
P,/mc Net 


11 图 (a): PEPI ARIA HEP 3E EI TE (p 1) f (p). 虚线 : 初始 分 布 , 点 ( 实 ) 线 : 采用 周期 性 (吸收 ) 边 界 条 件 的 PIC 模拟 获得 的 电 
子 在 Qet = 1200 时 刻 的 分 布 . 摘自 文献 [80]. 图 (b): 采用 不 均匀 ( 实 线 )、 均 匀 ( 虚 线 ) 背 景 磁场 的 电子 回旋 脉 泽 辐射 PIC 模拟 获得 的 被 激发 电磁 波 的 
能 量 随 时 间 的 演化 . 摘自 文献 [81]. 


E 


Fig.11 Panel (a): the electron perpendicular momentum distribution function f(p,) in different PIC simulations. The dashed 
line gives the initial distribution. The dotted (solid) line is from a PIC simulation with periodic (absorbing) boundary condition 
at Qet = 1200. See Ref. [80] for the original plot of panel (a). Panel (b): time history of the electromagnetic wave energy from 
the PIC simulations for the ECM emission with nonuniform (solid line) or uniform magnetic field (dashed line). See Ref. [81] for 


the original plots. 
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64 45 周 晓 伟 等 : 


Case B: y=1.2 (100keV), 9u=30 


空间 等 离子 体 中 相 3 


(d) 


=, = 


图 12 (a)-(d) 列 : 不 同 电磁 场 分 量 在 沿 95 方 向 的 


色散 空间 中 的 能 量 分 布 . 每 


(i) 列 展示 了 环 - 束 电子 能 量 的 减少 量 ( 一 人 x) 以 及 不 同 波 模 所 携 守 


带 的 能 量 随时 间 


F 射电 辐射 过 程 的 粒子 模拟 研究 


个 图 


3 期 


1000 1500 2000 


200 400 600 800 1000 1200 

tieġ] 
中 V 型 的 实 线 对 应 光波 . 每 个 图 的 颜色 代码 以 及 归 一 化 系数 都 相同 . 
的 演化 . 第 1、2 行 则 对 应 两 个 不 同 的 模拟 事件 . 该 图 摘自 文献 [84]. 


-80 -60 -40 -20 


Fig.12 Columns (a) to (d): dispersion spectra of different electromagnetic components along direction 05. The slope of V-like 


lines show the light wave. All panels have the same contour code and normalization. Column (i): time history of decrement of 


total kinetic energy (A Ex) of the ring-beam electrons, and energy histories of different wave modes. The first and second rows 


are for two different simulations. See Ref. [84] for the original plots. 


在 近来 的 系列 工作 中 ,Ning 等 B5-89 和 Youse- 
fzadeh 等 e] 的 PIC 模拟 研究 也 关注 强 磁化 (频率 比 


达 0.06% 和 0.17% ( 见 图 14). 而 损失 锥 电子 驱动 的 
回旋 脉 泽 过 程 主 要 激发 Z 模 和 X1 模 , X2 模 也 有 了 明显 


Re = 10 和 4) 等 离子 体 中 高 能 电子 的 激发 自由 
能 在 不 同 波 模 上 的 分 配 问 题 . 不 过 , 他 们 考虑 
到 X 模 基 波 在 传播 过 程 中 可 能 被 二 次 谐 波 吸收 


的 效应 , 在 其 研究 中 他 们 主要 关注 了 X 模 的 二 


谐 频 (X2) 的 激发 . 图 13、14 和 15 展 示 了 该 一 
拟 研究 的 主要 成 果 . 

BE, = = [E,, Ey, BRE, 
Bo 人 @) 方 向 上 的 色散 空间 的 


WB, = 


N 


yt 


ite (= 


系列 模 


中 ,图 13 通 过 不 同 电场 分 
ks /ks. 注 : 


ab EL 
H AR 
AER 


能 量 分 布 来 


展示 等 离子 体 中 不 同 波 模 的 激发 情 


了 不 同 电子 回 


旋 脉 泽 辐射 PIC 数值 


况 ; 图 14 则 展示 


模拟 获得 的 等 


增长 . 他 们 将 损失 锥 电子 驱动 的 X2 模 的 增长 归 因 
于 两 个 反 向 传播 Z 模 波 ( 即 Z + Z 一 X2) 或 者 Z 模 波 
与 X1 模 波 ( 即 Z + X1 一 X2) 的 非 线性 波 - 波 耦 合 过 
程 ( 见 图 15), 该 过 程 为 强 磁化 条 件 下 谐 频 辐射 产生 
的 新 机 理 , 与 Ni 等 "提出 的 Z 模 波导 致 激发 自由 
能 重新 分 配 的 想法 类 似 . 


5 ” 强 、 弱 磁化 等 离子 体 中 相干 射电 
辐射 机 制 PIC 模 拟 研究 的 比较 分 


离子 体 中 被 激发 波 模 的 能 量 随 时 间 的 演化 (图 中 各 
个 能 量 曲线 都 由 相应 的 高 能 电子 的 初始 动能 已 。o 


进行 归 一 化 ); 图 15 mr ed 
BUS Ak do = 及 来 展示 X2 模 生成 的 可 能 


dh 


对 比 前 面 关于 强 磁化 和 弱 磁 化 等 离子 体 环 境 
下 相干 射电 辐射 机 制 的 各 种 PIC 模 拟 结果 , 可 以 发 
现在 强 磁 化 等 离子 体 环境 下 高 能 电子 束 激发 自由 


IR] H6 


能 向 辐射 电磁 波 转换 的 效率 要 远 远 高 于 弱 磁 化 等 


原因 ,其 中 如 和 分 别 为 平行 TETERE 
(Bo = Boe) 方 向 上 的 波 数 . 他 们 的 模拟 结果 显示 : 


蹄 状 分 布 电子 和 损失 锥 电子 能 够 


效 激发 X2 模 的 


离子 体 环 境 中 的 转换 效率 . 这 说 明 频 率 比 Re > 1 的 


$ 


Hj, HEX2Ei JEDE BRA MBO AS PRI 


强 磁 化 等 离子 体 环 境 更 加 有 利于 相干 射电 辐 


E 5 E 


H 


不 同 的 机 种 


| 产生 . 他 们 指出 : 蹄 状 电子 可 通过 回 


磁 波 的 激发 和 增长 . 这 在 物理 


E 上 不 难 天 


Eit, KAE 


旋 脉 泽 不 稳定 性 直接 激发 Z 模 和 X2 模 ( 见 图 13)， 而 
且 X2 模 的 激发 效率 随 高 能 电子 束 对 背景 电子 的 密 


度 之 比 的 增 大 而 增强 , 对 于 np /no = 590110968] 
密度 比 , 高 能 电子 向 X2 模 的 能 量 转换 效率 分 别 高 


He 里 . 


34-15 


弱 磁 化 等 离子 体 环 境 中 ECME 机 制 不 能 有 效 工作 ， 
因而 激发 自由 能 没有 直接 转换 为 辐射 电磁 波 的 渠 
道 , 只 能 通过 非 线性 波 - 波 耦合 这 类 激发 自由 能 二 


; ^ Bed 
次 分 配 过 程 间 接 转 换 进 辐射 电磁 波 , 如 等 离子 体 辐 
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A: ATL tll P BA] BE ZR ZR BR PS T E A Jc I8] BB ZI E 
耦合 模式 . 而 强 磁化 等 离子 体 环境 为 BCME 机 制 有 
效 工作 提供 了 必要 条 件 , 高 能 电子 束 的 激发 自由 能 
可 以 通过 电子 回旋 脉 泽 不 稳定 性 直接 转换 进 辐射 
电磁 波 , 从 而 为 激发 自由 能 向 辐射 电磁 波 的 转换 提 
供 了 直接 的 转换 渠道 , 因而 辐射 电磁 波 能 够 获得 与 
静电 波 “ 平 等 "的 增长 条 件 和 能 量 分 配 机 会 . 

当然 , 等 离子 体 线性 不 稳定 性 的 理论 分 析 只 能 
提供 一 个 基本 的 参考 , 尤其 是 在 磁化 等 离子 体 的 


2 -10 -5 


复杂 环境 中 , 除了 静电 波 模 和 横 电 磁 波 模 之 外 , 还 
有 各 种 混杂 波 模 . 在 初始 的 线性 阶段 , 它们 有 可 能 
被 同时 激发 , 进入 非 线性 阶段 , 各 种 波 - 波 耦 合 过 程 
相互 竞争 也 可 能 影响 到 最 终 取 得 优势 的 饱和 波 模 . 
PIC 方 法 能 为 我 们 判断 最 终 的 优势 波 模 提 供 一 个 很 
好 的 参考 依据 . 不 过 , 前 面 提 及 的 模拟 研究 一 般 都 
是 选择 在 频率 比 R > 1 或 Re < 1 的 极端 磁化 环境 
中 进行 的 , 在 Re > 1 和 Re < 1 的 中 间 过 渡 磁 化 环境 
中 , PIC 模拟 研究 的 结果 显然 是 值得 期 待 的 . 
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13 HJ EE EI IRL ERO Ae MEE. 图 (a)、(b)、(c) 中 9 = 90°. AO (c))9 Jn Y Z KEREM) Of 
(c)) 的 理论 色散 曲线 (白色 点 线 ). 图 (a)、(b)、(c) 摘 自 文献 [85]. 图 (d) 和 (e) 展 示 了 电子 回旋 频率 (Q。。 = Wee) hi mu i E, TE IRI LE 

的 能 量 分 布 , 图 (全 展示 电场 分 量 互 , 的 能 量 在 频率 29.。< w < 2.259。。 范 围 的 分 布 . 图 (d)、(e)、(f) 摘 自 文献 [86]. 
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Fig.13 Wave dispersion diagrams of the electric components E; along direction 0. Panels (a), (b), (c) are with 0 = 90°. 
Analytical dispersion relations of the X/Z mode and O mode are overplotted as dotted white lines in panels (b) and (c), 
respectively. See Ref. [85] for the original plots of panels (a), (b), (c). Panels (d) and (e) display the diagrams of E, around 
Qee = Wee along different directions. Panel (f) displays the diagrams of E, around 2Q,. < w < 2.25€... See Ref. [86] for the 
original plots of panels (d), (e), (f). 
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图 14 不 同 的 电子 回旋 脉 泽 辐射 PIC 数 值 模拟 所 获得 的 被 激发 的 波 的 能 量 随时 间 的 演化 (图 (b) 摘 自 文献 [85], 图 (f) 摘 自 文献 [86], 图 (i) 摘 自 文 
献 [87]). 


Fig.14 Temporal variations of energies of the excited wave modes from different PIC simulations of ECME (panel (b) from Ref. 
[85], panel (f) from Ref. [86], and panel (i) from Ref. [87]). 
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图 15 图 (f): B5 EE,RIBERETER, — ki 空间 的 分 布 . 图 (a)-(e) 与 图 (f) 相 同 , (ELI COURTS C T- DX A E ME. 颜色 相同 的 箭头 用 于 展示 可 能 


的 Z + Z X2 Z4 XL > XQ. 黑色 箭头 则 对 应 于 所 生成 的 波 (X2). 该 图 摘自 文献 [86]. 


Fig.15 Panel (f): the energy distribution of E, in the kj — kı space. Panels (a) to (e) are the same with panel (f), but blue and 


brown represent results of E, and E, components, respectively. Possible wave-wave coalescence process Z + Z — X2, 


Z + X1 — X2 are overplotted with the same-color arrows. Black arrows represent the target mode (X2) to be generated. See Ref. 


[86] for the original plots. 


Zhou 等 B83 针 对 Re = 0.2. 0.3. 0.5. 1.0. 2.0. 
3.0 等 从 弱 磁 化 到 强 磁化 的 不 同等 离子 体 环境 中 各 
种 波 模 的 激发 情况 进行 了 PIC 模拟 研究 . 在 该 研究 
中 他 们 也 采用 了 环 - 束 型 动量 分 布 的 高 能 电子 束 ， 
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电子 束 的 相对 密度 比 为 wu/no = 5%. 图 16 展 示 了 
他 们 所 模拟 的 等 离子 体 中 被 激发 电磁 波 的 电场 分 
量 , 在 沿 0 的 色散 空间 k-w 空 间 的 能 量 分 布 (及, 为 
横向 电场 , 其 可 展示 等 离子 体 中 电磁 波 的 特性 ), 其 


64 卷 天 xX 


中 (a) 行 和 (b) 行 分 别 对 应 于 传播 角 9 = 20° ( 准 平行 
传播 ) 和 90* (垂直 传播 , 相对 于 背景 均匀 磁场 轧 = 
Boh) 的 情况 . 从 图 16 中 不 难 发 现 , 辐射 电磁 波 强 度 
随 频率 比 Re 增 加 而 增强 的 明显 趋势 , 在 Re < 1 的 情 
形 下 辐射 电磁 波 的 能 量 基本 可 忽略 不 计 . 这 一 结果 


学 报 3 期 


太阳 光斑 观测 中 , 相 较 于 硬 X 射 线 暴 , II 型 射电 暴 
一 般 会 滞后 < 1s 出 现时 9. 通常 倾向 认为 产生 太阳 
硬 X 射 线 暴 和 III 型 射电 暴 的 高 能 电子 是 源 自 同 

个 加 速 源 区 的 同一 批 电 子 束 , 只 是 由 高 能 电子 非 相 
干 轨 致 辐射 过 程 产生 的 硬 X 射 线 姑 一 般 来 自 相 对 致 


进一步 印证 了 前 面 分 别 在 强 磁 化 和 弱 磁 化 等 离子 
体 环 境 中 的 模拟 研究 结果 , 说 明 等 离子 体 磁化 程度 
在 相干 射电 辐射 机 制 中 确实 起 关键 作用 . 这 一 结论 
也 与 太阳 耀 斑 暴发 现象 中 经 常 观测 到 便 X 射 线 暴 、 
II 型 射电 暴 之 间 存 在 时 间 延 迟 的 现象 相 一 致 . 在 


WeelWpe = 0.2 WeelWpe = 0.3 WeelWpe = 0.5 


密 的 源 区 附近 , 一 般 为 弱 磁 化 等 离子 体 环 境 , 而 相 
干 射电 辐射 过 程 导致 的 II 型 射电 暴 则 需要 等 到 高 
能 电子 束 进 入 密度 较 低 的 强 磁 化 等 离子 体 区 域 才 
能 够 产生 . 


WeelWpe = 1.0 WcelWpe = 2.0 Wee/Wpe = 3.0 


9) 210910(|Ez/Bnorml) at 8 = 20.0° 


0077 15-10-05 00 05 10 15 00 15-10-0500 05 10 L5 00 15-10-0500 05 10 15 


0077 15-10-05 0.0 0.5 10 15 


» -15-1.0-0.5 0.0 05 10 1.5 n -15-10-05 0.0 05 10 15 


() 2iogao(|Ez/Bnom|) at 8 = 90.0 


UU LS-L0-05 00 0.5 LO L5 00 -15-10-05 00 05 10 L5 00 -15-10-05 00 0.5 L0 15 00 -15-10-05 00 05 L0 L5 09 -15-10-05 0.0 05 L0 15 % -L5-10-05 0.0 05 L0 15 
k [Wrorm / C] k [Wnorm | C] k [Wnorm | C] k [Wnorm / C] 


-8 7 -6 -5 EI 


3 E -1 0 1 


图 16 不 同 磁化 强度 (从 左 列 到 右 列 依次 为 Re = 0.2, 0.3, 0.5, 1.0, 2.0, 3.0) 的 等 离子 体 中 , WEEER ((a) 行 : 0 = 20°, (b)fr: 
9 = 90?) 的 色散 空间 的 能 量 分 布 . 每 幅 图 中 从 低频 至 高 频 受 加 展示 了 冷 等 离子 体 中 哨 声 模 ( 黑 色 点 线 )、Z 模 (黑色 虚线 )、O 模 (洋红 色 点 线 ) 以 


及 X 模 (洋红 色 实 线 ) 的 色散 曲线 . 图 中 Bnorm = Bo 和 wnorm = 5wpe 分 别 为 电磁 场 和 频率 的 归 一 化 常数 . 该 图 摘自 文献 [88]. 


Fig.16 Spectra of electric component E, along propagating direction 0 (row (a): 0 = 20°, row (b): 0 = 90?) in the dispersion 
space for differently magnetized plasmas (from left to right columns Re = 0.2, 0.3, 0.5, 1.0, 2.0, and 3.0). In each panel from low 


to high frequencies, the overplotted lines are the dispersion curves of the whistler (black dotted lines), Z (black dashed lines), O 


(magenta dotted lines), and X (magenta solid lines) modes in magnetized cold plasmas, respectively. Bnorm = Bo and 


Wnorm = 5Wpe in each panels are the normalization of electromagnetic field and frequency, respectively. See Ref. [88] for the 


original plots. 


众所周知 , 在 等 离子 体 线性 不 稳定 性 理论 中 ， 


的 高 能 电子 束 (nby/no = 5%) 在 磁化 等 离子 体 中 对 


除了 背景 等 离子 体 环境 的 磁化 状态 , 高 能 电子 束 
本 身 的 能 量 和 动量 分 布 也 能 够 对 各 种 波 模 的 激发 
效率 具有 重要 、 甚 至 关键 性 的 影响 . 考虑 到 太阳 和 
空间 高 能 电子 束 普遍 呈现 为 窜 律 能 谱 分 布 的 观测 
事实 80 9，Zhou 等 8B4 -95 采用 2.5 维 PIC 数值 模拟 研 
究 了 具有 各 向 异性 速度 分 布 和 低能 截止 究 律 能 谱 


辐射 电磁 波 的 激发 效率 以 及 其 相应 的 非 线性 饱和 
水 平 . 图 17、18、19 展 示 了 该 系列 模拟 研究 的 主要 
BR, FLU. Os ay cA sd WA Te BE Rt HL 
E^ {EG RE A LE EE TE. (AR BE ER) BE 
化 指数 、( 高 能 段 ) 谱 指数 、 投 射 角 中 心 和 投射 角 展 
宽 . 图 17 展 示 了 具有 低能 截止 容 律 能 谱 的 高 能 束 流 
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FA BAKU, ~ 0.5c 的 电子 束 激发 波 的 模式 以 静电 


IE x k/|k|). AlE = |B k/k) E W e E 
随时 间 的 演化 ( 注 : 图 17 (al) 和 (a2) 分 另 


I 展示 了 横 
向 电场 不 具有 和 具有 快速 增长 相 的 事件 ; 而 图 17 
(a3) 中 只 展示 了 具有 较 强 的 纵向 电场 激发 的 事 
fr, 即 事件 12 和 15). 图 18 展 示 了 图 17 中 事件 8 和 11 
的 横向 电场 | 蕊 | 在 沿 9 = 60* 的 色散 空间 的 能 量 分 
布 (第 1、2 列 , FER Po > 0 和 w < 0 的 区 域 分 别 
j 于 展示 了 右 、 左 旋 偏 振 的 电磁 波 ) 以 及 左 、 右 旋 
可 逃逸 电磁 波 的 能 谱 P(w)( 第 3 列 , 这 里 我 们 简单 将 
满足 w > wpe 和 |ck/w| < 1 的 电磁 波 定义 为 可 逃逸 
电磁 波 ). 图 19 则 通过 横向 电场 | 成 | 在 色散 空间 中 的 
能 量 分 布 来 展示 了 具有 低能 截止 守 律 能 谱 的 高 能 
束 流 电 子 对 不 同 磁化 强度 的 等 离子 体 中 电磁 波 的 
激发 的 差别 . 这 些 研究 结果 显示 : 在 Re > 1 的 强 磁 
化 等 离子 体 环境 中 , 寡人 律 谱 高 能 电子 束 的 低能 截止 
动量 必 能 够 显著 影响 静电 波 和 电磁 波 的 激发 优势 


— 


波 为 主导 , 而 当 0.1c < Uy < 0.5c 时 , 其 激发 的 电磁 
波 较 强 ( 见 图 17 (al) (a2) (a3)). 同时 , 高 能 电子 
束 的 容 律 能 谱 指 数 a 则 会 影响 被 激发 电磁 波 模 的 增 
长 率 和 饱和 强度 (图 17 (al)、(a2)、(a3))、 辐 射 传 
播 方 向 ( 见 图 17 (a)、(b)) 以 及 不 同 波 模 间 的 能 量 分 
配 ( 见 图 18). 另外 , 高 能 电子 束 速度 分 布 的 各 向 异 
性 , 尤其 是 投射 角 展 宽 1 对 辐射 电磁 波 的 激发 效率 
有 具有 明显 影响 : 投射 角 展 宽 越 窗 , 激发 辐射 电磁 波 
的 效率 越 高 ( 见 图 17 (al). (a2). (a3)). 针对 不 同 
频率 比 R = 2.5、1.0、0.5 的 模拟 研究 结果 则 进 一 
步 证 实 了 前 面 的 结论 : 在 ECME 机 制 起 主导 作用 的 
强 磁化 (Re > 1, 见 图 19 第 1 列 ) 等 离子 体 中 , 高 能 电 
子 束 所 激发 的 辐射 电磁 波 的 能 量 远 远 高 于 等 离子 
体 辐射 机 制 起 主导 作用 的 弱 磁 化 (Re < 1, 见 图 19 
第 3 列 ) 等 离子 体 中 的 辐射 电磁 波 能 量 . 
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能 截止 容 律 能 谱 的 高 能 束 流 电子 所 激发 波 的 横 ( 图 (al) 和 (a2))、 纵 (图 (a3)) 向 电场 的 能 量 随时 间 的 演化 . 其 中 
的 能 量 被 相应 事件 中 高 能 束 流 电子 的 初始 动能 skinetic-eeo 进 行 归 一 化 , 且 不 同 颜 色 的 曲线 用 于 区 分 不 同 的 事件 . 图 (a)、(b): 分 别 展示 了 事 


件 8 和 11 中 被 激发 的 沿 不 同方 向 传播 的 横向 电场 的 能 量 随时 间 的 演化 . 磁场 的 初始 能 量 =aorm 被 用 于 其 能 量 的 归 一 化 . 该 图 摘自 文献 [94]. 


Fig.17 Panels (al), (a2), (a3): energy evolution of the transverse (panels (al) and (a2)) and the longitudinal (panel (a3)) 


electric fields of waves excited by energetic cutoff power-law beam electrons, where energies are normalized by the initial kinetic 


energy of the energetic cutoff power-law electrons skinetic-ee0 in corresponding case. Different colors in each panels are used to 


distinguish different cases. Panels (a), (b): energy evolution of transverse electric field propagating along different angle bins in 


Cases 8 and 11, respectively, which are normalized by the initial magnetic field energy Enorm. See Ref. [94] for the original plots. 
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图 18 (1) (2)7l: 3 ES RU LHP HUE EZ | 在 沿 0 = 60° 的 色散 空间 的 能 量 分 布 . 不 同 的 列 对 应 不 同 的 时 段 . 在 w > 0 的 区 域 , 由 低频 到 高 频 分 

别 县 加 展示 了 冷 等 离子 体 中 哨 声 模 ( 灰 色 点 线 )、2Z 模 (灰色 虚线 ) 以 及 X 模 (灰色 虚 - 点 -点 线 ) 的 色散 曲线 . O 模 (灰色 虚 - 点 线 ) 的 色散 曲线 则 展示 
Fw < 0 的 区 域 . 所 有 的 图 共享 同一 颜色 代码 . (3) 列 : ERTA Ate n] 6356 ru UR DA) fei P (wo) , 在 不 同 w 区 域 中 不 同 颜色 的 实 线 用 于 区 分 
E(k > 0). (Kk < 0) 向 传播 的 电磁 波 . 为 了 对 比 同 向 传播 的 左 、 右 旋 可 逃逸 电磁 波 的 能 量 的 大 小 , 右 旋 可 逃 揭 电磁 波 的 能 谱 已 (w > 0)( 被 镜像 ) 也 
Eo Tw < 0 的 区 域 ( 见 绿色 、 黑 色 的 虚线 ). 该 图 摘自 文献 [94]. 


E 


Fig. 18 Columns (1), (2): energy distribution of transverse electric field E in dispersion space along direction 0 = 60° of 
Cases 8 (column 1) and 11 (column 2) over different time periods. In the w > 0 regions, from low to high frequency, the 
overplotted lines are the dispersion curve of the whistler (gray dotted line), Z (gray dashed line), and X (gray dash-dot-dotted 
line) modes in the magnetized cold-plasma limit, respectively. The dispersion curve of the O mode (gray dash-dotted line) is 
present in the w < 0 regions. The color-bar codes of each panel are the same. Column (3): power spectrum P(w) of the escaping 
right- and left-handed polarized waves. Different colored solid lines in each w region are used to distinguish the forward and 
backward propagating waves with k > 0 and k < 0, respectively. P(w > 0) of the escaping right-handed polarized waves are 
mirrored to the w < 0 region, shown by the green and black dashed lines, for a direct comparison for the energies of the escaping 


right- and left-handed polarized waves propagating along the same direction. See Ref. [94] for the original plots. 
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图 19 在 不 同 磁化 强度 的 等 离子 体 (Re = 2.5, 1, 0.5) 中 , WORCÉ IR HE A E REG HT e cR T PCR UII Eu LE? EAR FRE 


向 (9 = 50°,40° 以 及 30°) 的 色散 空间 中 的 


fitit (e ) 分 布 . 不 同 的 行 对 应 不 同 的 时 段 . 每 幅 图 中 全 加 展示 的 色散 曲线 与 图 18 相 同 . 另外 , 所 有 的 图 


Es [8] — t 


色 代码 . 该 图 摘自 文献 [95]. 


=> 
Fig.19 Energy distribution of transverse electric field |E;| (i.e., £g?) excited by energetic cutoff power-law electrons, in 


dispersion space along directions 0 = 50?,40? and 30? in differently magnetized plasmas (Rẹ = 2.5, 1, 0.5 shown in different 


columns), respectively, over different time periods (shown in different rows). The overplotted lines in each panel are the same 


with those in Fig.18. Additionally, the color-bar codes in each panels are the same. See Ref. [95] for the original plots. 
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6 总结 与 讨论 拟 研究 的 结果 之 间 也 存在 一 些 差异 , 这 些 差异 可 能 
正如 Treumann09 所 指出 的 , 射电 辐射 机 制 万 ee ado ide 
其 是 相干 射电 辐射 机 制 , 是 天 体 物理 中 最 复杂 的 电 tt a o: 
注 辆 射 机 制 ,一 方面 是 由 于 相干 辑 射 机 制 的 核心 过 |。 a ge mae hate 4 tn RM sh Te 
程 一 富 离 了 体 线性 或 /和 非 线性 不 稳定 性 的 激发 。 同 都 有 可 能 导致 得 到 明显 不 同 的 PIC 模拟 结果 
效率 十 分 敏感 地 依赖 于 源 区 等 离子 体 的 各 种 物理 Boos cu 
参量 TA AL, ER TSE BOG TABOR TR 发， 但 等 离子 体 集体 相互 作用 的 本 质 特 征 使 得 背景 
区 这 样 极 少 数 卫星 能 够 直接 进行 实地 探测 的 辐射 no es D A OLEH IE 
源 区 外 , 绝 大 多 数 空间 天 体 射电 辐射 源 区 的 物理 参 。 其 非 线性 耦合 的 微观 动力 学 过 程 , 这 就 要 求 PIC 模 
lL 只 能 通过 天 文明 测 的 手段 进行 间接 诊断 , 相应 的 ee a oo gues HY DICED 
诊断 精度 远 远 达 不 到 足以 唯一 确定 辐射 机 制 所 需 


, 由 于 离子 和 电子 在 微观 动力 学 特征 尺度 方面 存在 

要 的 精准 程度 , 因而 给 编 射 机 制 理论 模型 的 研究 带 gs aig pase JH I HHEUEOEPICHE IE 

”来 了 很 大 的 不 确定 性 . 另 一 方面 , 相干 射电 辐射 机 grew scien m Jhesu cR E 
pm 制 的 理论 研究 目前 只 是 在 线性 理论 的 研究 上 相对 内 的 风险 ， 所 以 一 个 成 功 的 PIC 模 拟 研 究 需要 在 


LO 比较 成 熟 , 线性 理论 研究 虽然 能 提供 许多 有 价值 的 模拟 计算 精度 和 有 限 计算 资源 之 间 寻 求 一 个 适当 


YT 启示 ， 但 我 们 实际 观测 到 的 相干 射 电信 和 号 还 需 经 过 的 平衡 . 从 现 有 的 PIC 模拟 研究 工作 来 看 也 确实 


c 非 线性 演化 并 达到 饱和 状态 的 结果 . 由 于 非 线 性 问 存在 一 些 需要 进一步 改进 的 空间 : 例如 ， 模拟 区 域 


© 题 在 数学 处 里 上 仍然 存在 许多 困难 ， 使 得 射电 辐射 所 采用 的 边界 条 件 . 有 关 的 模拟 研究 已 经 证 明 ， 不 
N 波 在 激发 增长 和 传播 阶段 所 遭遇 的 各 种 波 模 非 线 同 的 边界 条 件 会 导致 不 同 的 运算 结果 [eol, 如 何 设计 
性 而 合 方面 的 许多 问题 都 还 很 不 明朗 , 有 竺 进一步 ”一 个 合适 的 边界 条 件 , 是 进行 PIC 模拟 研究 必须 认 
澄清 . Wu HL, 实际 观测 到 的 天 体 相干 射电 辐射 都 已 真 考 虑 的 重要 问题 . 另 一 个 值得 注意 的 问题 是 高 能 
远离 了 其 辐射 源 区 , 与 线性 理论 预期 的 物理 特征 可 BE 子 束 所 携带 的 初始 电流 的 中 性 化 问题 3, 如果 不 
能 存在 显著 差异 , 这 也 给 天 体 相干 射电 辐射 机 制 的 适当 处 理 , 模拟 空间 中 将 出 现 一 个 强 的 净 电 流 , 该 
物理 诊断 增加 了 更 大 的 难度 . E 流 有 可 能 引起 某 些 波 模 , 尤其 是 静电 波 模 的 非 物 
针对 相干 射电 辐射 机 制 研究 遭遇 的 上 述 非 线 AK, 甚至 导致 模拟 结果 严重 失真 . 此 外 , 现 有 
性 理论 困难 , 基于 基本 的 粒子 动力 学 过 程 的 PIC 模 。 ”的 PIC 模 拟 研究 主要 是 对 相干 射电 辐射 源 区 有 关 物 
拟 方法 为 相干 射电 辐射 的 理论 研究 提供 了 一 个 很 。”” 理 过 程 的 再 现 , 而 实际 观测 到 的 射电 辐射 则 是 在 经 
— o0 好 的 研究 手段 . 本 文 介绍 了 近年 来 已 经 开展 的 有 关 ”过 了 一 系列 传播 动力 学 过 程 后 才 到 达观 测 者 的 , 在 
= 相干 射电 辐射 机 制 的 PIC 模拟 研究 工作 及 其 主要 结 这 个 漫长 的 传播 过 程 中 电磁 波 可 能 要 经 历 诸如 反 
R. 一 个 比较 一 致 的 结论 是 高 能 电子 束 在 Re > 1 的 射 、 折 射 、 吸 收 、 法 拉 第 旋转 等 诸多 物理 过 程 的 
强 磁化 等 离子 体 环境 中 激发 产生 的 辐射 电磁 波 强 。 ”扰动 和 转变 . 因此 , 为 了 能 够 将 实际 观测 结果 与 理 
度 远 远 高 于 在 Re < 1 的 弱 磁 化 等 离子 体 中 被 激发 。 ” 论 预期 进行 更 好 的 对 比 、 诊 断 , 我 们 也 需要 对 相干 
的 辐射 电磁 波 的 强度 . 在 强 磁 化 等 离子 体 中 , 相干 射电 辐射 在 行星 际 或 星际 空间 的 传播 过 程 进行 相 
射电 辐射 机 制 由 ECME 机 制 主 时 , 而 在 弱 磁 化 等 离 ”应 的 PIC 模 拟 研究 986-93]. 

子 体 中 则 由 等 离子 体 辐射 机 制 主导 . 所 以 , PIC 模 目前 , 空间 等 离子 体 中 相干 射电 辐射 机 制 
拟 研 究 的 上 述 结论 也 意味 着 : 相 比 于 基于 非 线性 的 PIC 模拟 研究 还 处 在 一 个 初始 的 发 展 阶 段 ， 随 
波 - 波 耦 合 过 程 的 等 离子 体 辐 射 机 制 , 基于 等 离子 着 模拟 程序 的 不 断 改进 和 计算 能 力 的 持续 提升 , 我 
线性 不 稳定 性 激发 的 ECME 机 制 能 够 更 加 有 效 地 们 相信 这 一 研究 领域 会 有 更 多 、 更 好 的 研究 结果 
激发 等 离子 体 中 的 电磁 波 辐射 . 然而 , 不 同 PIC 模 值得 期 待 . 
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Probing Coherent Radio Radiation Process in Space Plasmas via 
Full-kinetic Particle-in-Cell (PIC) Simulations 


ZHOU Xiao-wei WU De-jin CHEN Ling 
(Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023) 


Asstract The radio emission mechanisms, especially the coherent radio emission mechanisms of radio 
bursts, are one of the most complicated and controversial topics in astrophysics. Due to the complex 
influence of multiple physical factors, theoretical studies of coherent radio emission mechanisms are very 
difficult. There has been a constant controversy between the two types of coherent emission mechanisms, 
i.e., the plasma emission and the electron cyclotron maser emission, for a long time. In recent years, people 
have been trying to apply the fully-kinetic Particle-in-Cell (PIC) numerical simulations to the study of 
coherent radio emissions, and have made some positive progress. In addition to reviewing some existing 
problems and difficulties in this topic, this paper will focus on the research of PIC simulations for the 
coherent radio emissions from space plasmas and the main progress in recent years. 


Key words plasmas, waves, radiation mechanisms: non-thermal, methods: numerical 
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